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Abstract 

We numerically study the vortex core structure in a noncentrosymmetric supercon- 
ductor such as CePtsSi without mirror symmetry about the xy plane. A single vortex 
along the z axis and a mixed singlet-triplet Cooper pairing model are considered. 
The spatial profiles of the pair potential, local density of states, supercurrent den- 
sity, and radially-textured magnetic moment density around the vortex are obtained 
in the clean limit on the basis of the quasiclassical theory of superconductivity. 

Key words: CePtsSi, Unconventional superconductivity, Vortex core, Local 
density of states, Core magnetization, Broken inversion symmetry 
PACS: 74.20.Rp; 74.70.Tx; 74.25.Op 



Much attention has been focused on the heavy fermion superconductor CePtsSi, 
which has a noncentrosymmetric crystal structure without mirror symmetry 
about the xy plane [1]. CePtsSi is an extreme type-II superconductor and the 
vortex structure of the mixed state in this system was recently studied by 
Kaur et al. [2] and Yip [3] on the basis of the Ginzburg-Landau theory and 
the London theory. In this paper, we investigate the vortex core structure on 
the basis of the quasiclassical theory of superconductivity [4], which enables 
us to calculate more microscopically the physical quantities such as the pair 
potential, local density of states, supercurrent density, and magnetic moment 
density. We consider a single vortex along the z axis in the clean limit. 

The noncentrosymmetricity (or the lack of inversion symmetry) leads to the 
mixture of Cooper pairing channels of different parity [5]. We consider the 
following superconducting order parameter in a singlet-triplet mixing form: 
A k = (tyfro + d k ■ cr)i6~y = ty&o + A(— k y a x + k x a y ) ia y , with the s-wave 

pairing component $ and the d vector d k = A(—k y ,k x ,0). This s + p-wave 
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pairing state is proposed for CePtsSi in Ref. [6]. Here, (a x , a y , a z ) are the 
Pauli matrices in the spin space, <r is the unit matrix, and k = (k x , k y , k z ) = 
(cos (f) sin 9, sin sin 9, cos 9) . 

The lack of inversion symmetry here is incorporated through a Rashba-type 
spin-orbit coupling with a form proposed in Ref. [7]. It splits the Fermi surface 
into two ones (I and II) by lifting the spin degeneracy [6]. From the original 
Eilenberger equation for noncentrosymmetric superconductivity [8], we obtain 
two equations corresponding to these split Fermi surfaces I and II in the case 
of the above s + p-wave pairing state [9] , 

ivi,n ■ V<?i >n + [iuj n f 3 - Ai >n , <?i,ii] = 0, (1) 

where Ayi = [(fi+if2)Ai 5 n — (fi— if2)Aj n ]/2, Ayi = \P ± A sin 9 are the order 
parameters on the Fermi surfaces I and II, (fi, f 2 , f 3 ) are the Pauli matrices in 
the particle-hole space, and the commutator [a, b] = ab — ba. We neglect the 
vector potential in Eq. (1) assuming the extreme type-II superconductivity. 
We use units in which h — k-Q — 1. 

The Green functions gi y u on the Fermi surfaces I and II are written as a matrix 
in the particle-hole space, 

!te(r,k,iu> n ) = -i«( *f ifl n U ) 

V — * Ji.ii — 51,11 / 

The regular Green function g as a matrix in the spin space is given by [6,8,9] 

g = gi°i + gimi 
= \( _gi + gn -K(9i-gn)\ 

2 \ -k'_(gi-g u ) gi + gn ) ' 

with aiji = (o"o ± Qk ■ o")/2 and = (—k y ,k x ,0). Here, k'± = k y ± ik x and 
k = (k x , k y , 0) = (cos (/>, sin 0,0). 

We consider a single vortex which has a form, Ai t n(r, <p r ; 9) = [^(r)±A(r) sin 9} exp(i0 r ). 
Here, the real-space coordinates r = (r cos0 r , r sin r , 0), and the vortex cen- 
ter is situated at r = 0. The Fermi surface is assumed to be spherical, and the 
differences of the density of states and the Fermi velocity v^ u between the two 
Fermi surfaces I and II are assumed to be small and are ignored. The results in 
this paper depend predominantly on the spin structure [Eq. (3)] and the gap 
structure on the 3D Fermi surfaces, and the Fermi-surface anisotropy would 
not lead to qualitatively different results as long as the spin and gap topologies 
are not altered. We numerically solve the gap equations given in Ref. [6,8,9] 
and the Eilenberger equations in Eq. (1) self-consistently as in Ref. [10]. When 
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(3) 
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solving the gap equations, we adopt the same values of parameters as used in 
Ref. [8]. Thus, both the pair potentials A and ^ are real and positive, and 
|A| > 1^1 [8]. From now on, T c is the superconducting critical temperature 
and £ = v-p/T c is the coherence length at zero temperature (v F = \v F \ is the 
Fermi velocity). 
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Fig. 1. Spatial profiles of the pair potentials. T/T c = 0.1-0.9 from top to bottom 
by 0.1 step, (a) the p-wave component A, and (b) the s-wave one ^. 



In Fig. 1, we show the spatial profiles of the pair potentials A (p-wave com- 
ponent) and ^ (s-wave one) around the vortex for several temperatures T. It 
is noticed that while the amplitude is different between A and \I>, the charac- 
teristic recovery length (namely, the core radius) is the same for both. 




Fig. 2. The density of states in the bulk without vortices. T/T c = 0.1 and 
r) = O.OlTc. 

The local density of states (per spin) is calculated by 



N(E, r) = -^Re(Tr[g(iu n + ir,)]) 

= ^Re(g I + g u ), (4) 

where (• • •) denotes the average over the Fermi surface, N is the density of 
states per spin at the Fermi level, and rj (> 0) is the energy smearing factor. 
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Before going into the vortex bound states, let us see in Fig. 2 the density of 
states in the bulk without vortices. There are four gap edges (solid line). The 
system has two split Fermi surfaces I and II [8,9]. The two of the gap edges 
originate from the fully-gapped Fermi surface I (dashed line), and the other 
two originate from the line-node-gap Fermi surface II (dash-dotted line). 




Fig. 3. The local density of states N(E,r) inside the vortex core in units of No. 
T/T c = 0.1 and rj = 0.05T C . Large zero-bias peak at (E,r) = (0,0) is truncated in 
the figure. 

In Fig. 3, we show the local density of states inside the vortex core. There 
are four branches of peaks, which are related to the vortex bound states. 
The outer (inner) two branches originate from the vortex bound states of the 
quasiparticles on the Fermi surface I (II). Thus, the present spectra inside the 
vortex core in the clean limit possess the same structure as those in a two-gap 
superconductor. 
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Fig. 4. Spatial profiles of the supercurrent density \j\ in units of 2\6\vfNqT c . 
T/T c = 0.1-0.9 from top to bottom by 0.1 step. Inset: \j\ at T/T c = 0.1 (solid 
line), the contribution of the Fermi surface I (dashed line), and that of the Fermi 
surface II (dash-dotted line). 

In Fig. 4, we plot the supercurrent density which is calculated by 



j = eT'£ i No(v F Tr[& (-iTrg)] 
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= -iireT^2N (v F (g l + gn)), 



(5) 



where e is the electric charge of the quasiparticle. We have confirmed numer- 
ically that \j\ decays as ~ 1/r far away from the core. \j\ exhibits essentially 
the same structure as that in usual s-wave superconductors. 
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Fig. 5. The radial component of the magnetic moment density M around the vortex 
in units of ii b NqT c . T/T c = 0.1-0.9 from bottom to top by 0.1 step. Inset: M at 
T/T c = 0.1 (solid line), the contribution of the Fermi surface I (dashed line), and 
that of the Fermi surface II (dash-dotted line). The minus value means that the 
sense of M is toward the vortex center. 

Finally, we investigate the magnetic moment density M. The vortex-core mag- 
netization in the present noncentrosymmetric system has been reported by 
Kaur et al. [2] and Yip [3]. Here, we calculate it to obtain numeric results 
at various temperatures by means of a more microscopic derivation. M is 
calculated by 



M = l i B TY i Nom*(-mg)\), 



(6) 



where fj, B is the magnetic moment of the quasiparticle. Substituting Eq. (3) 
into this, we obtain 



M x = -iTc/j, B Tj2 N o((-k y )(9i ~ 9n)), 

M y = -iirpi B Tj2No(k x (gi ~ 9n)), 
M z = 0. 



(7) 

(8) 
(9) 



We have confirmed numerically that the azimuthal component of M is zero 
around the vortex, namely M is aligned in the radial direction in the xy 
plane. In Fig. 5, we show the spatial profiles of the radial component of M. 
We checked that \M\ decays as ~ 1/r far away from the core. The spatial 
profiles of \M\ are similar to those of but there is a difference because of 
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their microscopic origin. While \j\ is composed of g\ + gu in Eq. (5), \M\ is 
composed of g\ — gu in Eqs. (7) and (8) [see also the insets of Figs. 4 and 5]. 

In conclusion, we calculated the vortex core structure in the noncentrosym- 
metric superconductor. We found that the singlet and triplet components of 
the order parameter have the same characteristic length (core radius) as each 
other. The local density of states exhibits the two-gap property in the clean 
limit, which is expected to be observed by STM [11,12] in order to clarify 
the validity of the present mixed-parity pairing model. We mention here that 
impurity effects in the noncentrosymmetric system can be different from those 
in simple two-gap systems in general, therefore interesting phenomena in the 
density of states are expected when the system deviates from the clean limit. 
Detailed investigations on such impurity effects are left for future studies. The 
radially-textured magnetic moment appears around the vortex, which is ex- 
pected to be observed directly by spin-polarized STM [13] in order to clarify 
the electronic structure [Eq. (3)] specific to the noncentrosymmetric supercon- 
ductivity with the rotating spin on the Fermi surfaces as shown in a picture 
in Ref. [14]. 
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